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The strong and electromagnetic decays of X(4350) with quantum numbers J p = ++ and 
2 ++ have been studied by using the effective Lagrangian approach. The coupling constant 
between X(4350) and D*D* is determined with the help of the compositeness condition 
which means that X(4350) is a bound state of D*D* . Other coupling constants applied in 
the calculation are determined phenomenologically. Our numerical results show that, using 



the present data within the present model, the possibility that V(4350) is a D*D* molecule 
can not be ruled out. 
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I. INTRODUCTION 

Recently, a hidden charm resonance named A"(4350) was observed by Belle collaboration in the 
analysis of the 77 —> 4>J/ ip process [l|] . The mass and natural width of this resonance are measured 



> 

t> 

to be (4350.6±|-f (stat) ± 0.7(syst)) MeV/c 2 and (13.3^ 7 f 9 (stat) ± 4.1(syst)) MeV, respectively. 
\Q ■ 

The product of its two-photon decay width and branching fraction to (f>J/ifi is (6.7^2 "4 (stat) =t 

O ! 

l.l(syst)) eV for J PC = 0++, or (1.5lg;£(stat) ± 0.3(syst)) eV for J PC = 2++. In literature, the 
structure of X(4350) has been proposed to be cess teraquark state with J PC = 2 ++ M, D*D* n 
% ! ^ s tate a L P-w»e Stete A fl. And — g th e ^Jj^ 

of the final states J/ip<j>, X(4350) can also have quantum numbers J PC = 1 l ". In Ref. [5], it 
was shown that A"(4350) cannot be a 1 ^ exotic D*D* molecular state. In this paper, we will 
accept it as a bound state of D*D* to study its strong and electromagnetic decays in the effective 
Lagrangian approach in case of J PC = ++ and 2 ++ . 

Since the mass of X(4350) is about 80 MeV below the threshold of D*D* s0 (m^ = 2317.8 ± 
0.6 MeV and m t = 2112.3 ± 0.5 MeV \m), it is reasonable to regard X(4350) as a bound state of 
D* D*. And because the quantum numbers of D* and D* are J p = + and J p = 1~ respectively, 
to form a bound state with quantum numbers J PC = ++ or 2 ++ , the coupling between X(4350) 
and its constituents should be P— wave. To determine the effective coupling constant between 
Y(4350) and it constituents D*D* , as in our previous work (for example Ref. |7j), we resort to the 



compositeness condition Z v = 0(Z Y as the wave function renormalization constant of A(4350)) 

n q 

which was early used by nuclear physicists [8j, [9( and is being widely used by particle physicists (see 
the references in Ql ) . Recently, this method has been applied to study the properties of some 



"exotic" hadrons 7|, 



10H18J] and some conclusions were yielded comparing with data. For other 



interactions, we write down the general effective Lagrangian and determine the coupling constants 
with help of data, theoretical calculation, SU(4) relation or the vector meson dominance (VMD). 

nnn 

As in our previous work [7|, llOHlg] , we introduce a correlation function including a scale parame- 
ter A x to illustrate the distribution of the constituents in the bound state A (4350). The parameter 
A x is varied to find the physical region where the data can be understood. In the physical region 
of A x , the partial widthes for strong and electromagnetic decays are yielded. 

This paper is organized as the following: In section [II] we will provide the theoretical framework 
used in this paper. We will present the analytic forms for the radiative and strong decay matrix 
elements and partial widths of A(4350) in section HTT1 And, the last section is our numerical results 
and discussions. 



II. THEORETICAL FRAMEWORK 

In this section, we will propose the theoretical framework for the calculation of the strong and 
electromagnetic decays of A (4350). 

A. The Molecular Structure of X (4350) 

As was mentioned above, we regard A(4350) as a D*D* si) bound state. And concerning the 
experimental status, we accept the quantum numbers of X(4350) as J p = ++ and 2 ++ . For 
scalar case, one can write the free lagrangian of A (4350) as 

£free = ^Xd^X - X -m\X 2 , (1) 

with m x as the mass of A (4350). The propagator of A (4350) can be easily written as 

Gf(x) = [ - 2 -e~^, (2) 

J (27rj 4 p z — m x — le 

which satisfies 

(d 2 + m 2 )G F {x) = -i5^(x) . (3) 
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While for tensor resonance we have the free Lagrangian as [1 

£free = --jX^vD^^Xxa , (4) 

where the symmetric tensor X^ v = X Vil denotes the J PC = 2 ++ field for X (4350) and 

+g Xa d^d u + g^d x d a - l -(g va d^d x + g uX d^d a + g^d v d x + g^ x d v d a ) . (5) 

The propagator for X M! ,(4350) is obtained as 

n x (x) = f - P ^ e~ ip ' x 

J (zvrj 4 — mc. — ze 

P - _„ +Mfi 

D^ Xa G x ^ = -i-{g^ a g v ^ + g va g^)5^{x) . (6) 

With respect to the discussions given in first section, one can write the effective Lagrangian 
describing the interaction between X(4350) and D*D* as 

Cf nt = -^=g s X(x) J dx 1 dx 2 C fJifl (x 1 ,x 2 )$ x ((xi - x 2 ) 2 )5(x - uj v xi - u s x 2 ) , 

i f r l 

£j£ t = -j=g T X^ LU {x) J dxidx 2 C tiv (x 1 , x 2 ) + C^(xi, x 2 ) - jg^Caa(xi, x 2 ) 

x$x((xi - x 2 ) 2 )5(x - uj v x\ - lo s x 2 ) , (7) 

where £j^ t is for scalar resonance case while Cf nt is for tensor resonance case. g s and g T are the 
effective coupling constants for the interaction between X (4350) and D*D* in scalar and tensor 
resonance cases, respectively. oj v and w s are mass ratios which are defined as 

uv = r 2 — ' Ws = r 29 — • ( 8 ) 

m nt +m n , m n ,+m n , 

^x{{x\ — x 2 ) 2 ) is a correlation function which illustrates the distribution of the constituents in the 
bound state. Fourier transform of the correlation function reads 

<My 2 ) = / ^*x(p 2 )e-^ y . (9) 

To write down Lagrangian 0, for simplicity, we have defined the tensor C^ u as a function of the 
constituents with the explicit form 

C M ,(x 1 ,x 2 ) = D* s + ( Xl )d v D^{x 2 ) + D* s -{ Xl )d^Dl + {x 2 ) . (10) 



The coupling constants g s and g T can be determined with help of the compositeness condition 
Zx = with Zx as the wave function renormalization constant of X(4350) which is defined as the 
residual of X(4350) propagator, i.e., 

d 



Z X = l-g x ^x( P 



(11) 



where g x = g s for scalar case while g x = g T for tensor case. For scalar resonance X(4350), 
g 2 T, s (p 2 ) = H s (p 2 ) is its mass operator. But for tensor resonance X(4350), g 2 T lj ,(p 2 ) relates to its 



mass operator via relation 



n^V) = \(a*x9vp + g^9»*)g 2 T z T ( P 2 ) + • • • , 



(12) 



where " • • • " denotes terms do not contribute to the mass renormalization of X(4350). The mass 
operator of X(4350) is illustrated by Fig. [TJ 



D* 



X(4350) 




X(4350) 



FIG. 1: The mass operator for X(4350). 

Concerning the Feynman diagram depicted Fig. Q] one can calculate the mass operator explicitly. 
To get the numerical result of the coupling constant g x , an explicit form of <&x{p 2 ) is necessary. 
Throughout this paper, we take the Gaussian form 

<S>x(p 2 ) = exp(p 2 /A 2 x ) , (13) 

where the size parameter A x parametrizes the distribution of the constituents inside the molecule. 
In the following calculation, we will find the physical value of A x by comparing our calculation of 
the product of X(4350) to two-photon partial width and the branching fraction to J/ip<f>. It should 
be noted that choice (fl~3]) is not unique. In principle any choice of &x(p 2 ), as long as it renders the 
integral convergent sufficiently fast in the ultraviolet region, is reasonable. In this sense, ^x{p 2 ) 
can be regarded as a regulator which makes the ultraviolet divergent integral well defined. 

With these discussions, we can calculate the effective coupling constant g x numerically. In the 
typical nonperturbative region A x = 1.0 ~ 2.0 GeV, using the central value of X(4350) mass, our 
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numerical result is found to be 

g s = 31.49 ~ 15.19 , g T = 62.70 ~ 34.54 . (14) 

In Fig. [2] we plot the A x dependence of the coupling constants. One can see that both coupling 
constants decrease against A x . This can be understood from the momentum integral of the mass 
operator. For scalar X(4350), the loop integral is quadratically divergent so the derivative of 
the mass operator which proportional to the inverse of g 2 increases against A x which means the 
coupling constant g s decreases against Ax- Similar argument can be given for g T . 
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FIG. 2: A x dependence of the coupling constant g x . 

B. Effective Lagrangian for Strong and Electromagnetic Decays of X(4350) 

The effective Lagrangian for the study of strong and electromagnetic decays of X(4350) consists 
of two parts: the electromagnetic part C cm and the strong part £ s t T . 

The electromagnetic interaction Lagrangian C em includes five parts: £^ from the gauge of 
the nonlocal and derivative coupling of Eq. (|T|), £f^ ge from the gauge of the kinetic term of the 
charged constituents D* and D*, the electromagnetic interaction Lagrangian including D* 
and D*, for electromagnetic interaction including D sl and D* and for electromagnetic 
interaction including D s \ and D* Q . 

One can write £^ by substituting C^ v with C^ n&c in Eq. ([7|) with 

C^ C (x 1 ,x 2 ) = e- ie/ ^^^) J D: + (x 1 )(^ + ze^(x 2 ))^V(x 2 ) 
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-J PC =2" 
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+e ieI(x uX r,P) D *- {xi) ^ _ ieA ^ X2 ))D* s + (x 2 ) , (15) 

where the Wilson's line I(x, y, P) is defined as 

I(x,y;P) = f X dz^(z) . (16) 

In our following calculation, the nonlocal vertex with one-photon is necessary. The nonlocal vertex 
with one-photon comes from two sources: One is from covariant derivative and another one is from 
the expansion of the Wilson's line. One can easily derive the Feynman rule for the nonlocal vertex 
with one-photon which comes from the covariant derivative. But to derive the Feynman rule for 
photon from Wilson's line, one may use the path-independent prescription suggested in 



20, 



21|. 



The electromagnetic vertex Cem gc from the gauge of the kinetic terms of the charged constituents 
can be easily written as 

£fr ge = ieA»(D* s0 -d^D* s +) + ieA IM [-D* s rJ fi D* s .+ + \d*~XdZ + ^*^*+] . (17) 

One can generally write the effective Lagrangian for electromagnetic interaction including D* 
and D* as 

£fm = eg DtoD *jV- u D* s0 + - V+D* s0 -]F^ . (18) 

where V~, is the gauged field strength tensor for D*^ with definition V^ v = (d^ =F ieA^D*.^ — 
{d v =F ieA u )D*.^. And similarly, the general effective Lagrangian and can be written as 

£em = e 9D sl D*^^ ua p[D~ 1 .^Dl+-D+ 1 .^D*.~}F al 3, 

£el = - ie 9D sl D* s0 ~fe^ al 3[D* s + D~ v ^ - D* s0 D^.^Fap . (19) 

Similar as the definition of V^ ul we have defined the gauged field strength tensor for with 
definition = (d, T ieA,)D* sl % - (d u T ieA v )D* sl %. 

The relevant coupling constants can be determined phenomenologically. Confined by the ex- 
perimental status, one cannot fix q^, from data, so we turn to the theoretical calculations (for 
example Ref. [7] and references therein). From literature, one can see that the minimal result of 
the theoretical calculation of D* —¥ D*^ decay width is 0.2 KeV. From this decay width, we get 
o , n . > 3.02 x 10~ 2 GeV^ 1 . 

The coupling constants <7d s1 d* 7 an d 9D sl D* s0 -y can be determined by using the HQET and branch- 
ing ratio for relevant processes. First, consider the decay of D s \ — > D s j, the effective Lagrangian 
can be written as 

£?m li?s7 = ieg DslD ADtD^ ^ - DJD+. ^]F, U , (20) 



where D s \- ^ v = d^D s i ]U — d u D s i-^ and F^ u = d^A u — d u A^. From this Lagrangian, one can express 
the decay width as 

r (Dil ^ D , 7 ) - ^K.,-<f. ( 21 ) 

Dsi 

The numerical result of the decay width has been evaluated by several groups. From the references 
given in Ref. [u]], we see all the results are larger than 0.6 KeV. So that we have gD sl D 3 -y > 
2.67 x 10~ 2 GeV -1 . The coupling constant 5D si D*7 relates to gD 3l D a -y via HQET as 

gn^D^ = 1 = 3m x 1Q -l QeV -i (22) 

9D sl D* sl m Dsl jrrvBl 

so that we have <7d s1 d*7 = 6.81 x 1CT 2 . The coupling constant g „ can be determined by using 
the relevant branching ratio given in PDG [g]. From f)19|) we have 

F(D sl D* o7 ) = ^|§^(mL - m^) s • (23) 

Using the central value of the branching ratio we have T(D s i — > D* "f)/T(D s i — > D s j) ~ 0.21 
which leads to g _ n . = 3.53 x 10" 2 GeV" 1 . 

In addition to the Lagrangian (J7|), the strong part £ s tr involves VVV-type Lagrangian describing 
the interaction of three vector mesons, the SW-type Lagrangian describing the interaction of one 
scalar meson with two vector mesons, SSV-type Lagrangian describing the interaction of two scalar 
mesons with one vector meson, AVV-type Lagrangian for the interaction of axial-vector with two 
vector meosns and .AS V- type Lagrangian for axial- vector-scalar- vector meson interaction, i.e., 

+ig m [D* s .-(D* s A(f>v) + D* + (<j> v d^D*-) + <P,(D* s -d^D* s +)] , (24) 
4 5 t r = W D . ID; -D:+ U - D*+D* s r^ v + g [D: -D:+ u - D*+D* s r^ v , (25) 

4 5 tf y = D . MD* s o~Wa + ) ~ i9,n* D *M D *^ D *» +) ' (26) 

^ sO su 7 sO s(l 

-9^D i D sl W/3 i D 7l; ^Zjt ~ D tl; yPl^afi » (27) 



£str V = -ig^D* D sl e,i»\o[D* Q ip^D+.^-Dl+^D 



si: Act J 



-ig<t>D* D sl e^Xa [D* s0 ^uD+. Xa - D* s +^ v D sl . A J . (28) 

Because of our less knowledge, we can not determine these coupling constants from data. Here, 
we resort to the vector meson dominance (VMD) model j3] . In the VMD model, the virtual 
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photon in the process e D*^ — > e D*^ is coupled to vector mesons 4> and J/ip, which are then 
coupled to D*q~. For zero momentum transfer, one has relation 

E J = e ' (29) 

V=</),i/i V- 

where "y v is the photon-vector-meson mixing amplitude 

C'V— 7— mixing = Tv^M^/U ' (^0) 

which can be determined from V — > e + e~ decay width, i.e., 

Tvee - , (31) 

where we did not include electron mass since it is much smaller than vector meson mass. For eft 
meson, using F(<f) -»• e+e") = 2.97 x 1(T 4 x 4.26 MeV fl| we have 7^ = 23472.3 MeV 2 , while 
7^ = 259965.8 MeV 2 when -> e+e - ) = 5.94% x 93.2 KeV is applied. Concerning that the 
virtual photon sees the charge of charm quark in -D* + meson through ipD* D* coupling and the 
charge of anti-strange quark in D*^~ meson through <pD* D* coupling, we have relations 

7v>SW r n* 2 7*9 0D . D . 1 

y sO sO g v sO sQ g f32) 

m 2 3 m 2 3 

From these relations we have , n » „* = 7.45 and o » „ = 4.47. To determine coupling constants 
5,/n* n* > we make extension to the VMD model used above, i.e., substituting Eq. (|29l) with 

V SO s 

—^2 = > ( 33 ) 

V=0,-0 V 

where D{ and -Dj denote the relevant charmed-strange mesons. Similarly, Eqs. (|32ft should also be 
extended to 

l^DiDj 2 J^DiDj 1 

— -2 - oWDiDj-y , — —n -^eg DiDjl . (34) 

m V) 6 m 4> 6 

From which we yield the relevant coupling constants as 

tyDiDj = 7.45 x g D . D . n , g4>D . Dj = 4.47 x g D . D .^ . (35) 
To fix the magnitude of coupling constant gy £>•£>•, we resort to the 517(4) relation as was used 



in Ref. 



23J] from which we have relations 



2 

9iPD*D* = -^9<j>D*D* = g-ipD* D* = 7.64 . (36) 



To fix the relative sig ns for the relevant effective Lagrangian, one can use HHChPT including 
D* and D s \ mesons [24[]. But even with this consideration, the relative signs of and C^ v to 
the other terms cannot be determined. We leave this as an ambiguity and discuss different cases 
in the following calculation. In summary, our framework of the interaction Lagrangian is 

r — r NL _l rgaugc , r SV , r AV , r VVV , r SVV , r SSV , r AVV , r ASV (o 7 \ 

Cg* = ±[£g + £tri (38) 

Up to now, we have fixed all the coupling constants that are necessary for our following calcu- 
lation of the electromagnetic and strong decays of X(4350). 



III. ELECTROMAGNETIC AND STRONG DECAYS OF X(4350) 

In this section, we will present the general forms of the matrix elements and partial widths 
for the electromagnetic and strong decays of X(4350) and the Feynman diagrams included in our 
calculation. 



A. Electromagnetic decay of X(4350) 

The four kinds of diagrams depicted in Fig. [3] and their corresponding crossing ones should be 
taken into account to study X(4350) — > 2j decay. Diagrams (A) and (B) are from the final state 
interaction due to the exchange of D*, D s \ and D* , diagram (C) arises from the gauge of the 
nonlocal and derivative coupling between X(4350) and its constituents D*D* but diagram (D) is 
from the Lagrangian (|18j) . 




FIG. 3: Feynman diagrams for decay X(4350) — > 77(cross diagrams should be included). 

For X(4350) with quantum numbers J PC = ++ , concerning the C/(l) em giauge invariance and 
the transverseness of the photon polarization vector, one can write down the matrix element for 
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the decay of X — > 2 7 as 



iM| m 



9a/3 



qi ■ Q2 



ea(<?i)e/j(g2) 



(39) 



While for tensor meson X(4350) with quantum numbers J FC = 2 ++ , its polarization vector satisfies 
e^ v = e u ^ and e£J = 0, so that the matrix element for electromagnetic decay can be written as 251.126 1 



ie 2 j-F^j^ 



9o 



Q2aQl/3 



7 (2) 
T-J-27 



qi ■ <?2 



(40) 



)(9"P ~ + fart ~ J Z^-)i.9v 



q* q 
where q = q\ — q 2 - From Eqs. ()39|40|) we express the decay width for X(4350) as 

2tt 



^e^ u (p)e a {qi)ep{q 2 ) 



r s (x -> 2 7 ) 



m x 

7T 



rv 2 F 2 



T T (X -> 27) = ^«L(54T-.2 7 " 4i^ 27 i^U 27 + 32F^ 27 ) , (41) 
where the subindices "S" and "T" denote the scalar and tensor resonance X(4350), respectively. 



To get the last equation, we have applied the sum of the polarization vector for tensor meson 



27] 



(P) 

polar 



1 



1 



(42) 



where 



(IV 



-3^ + l^r- 



B. Strong decay of X(4350) 



We should take into account the Feynman diagrams illustrated in Fig. [J] in the study of the 
strong decay of X(4350) —> J/ip<f>. Furthermore, in addition to these diagrams, their crossing ones 
should also be included. 



X(4350) 



J/4' 
D*JD sl 

Vwv 



X(4350) 




(A) 



(B) 



FIG. 4: Feynman diagrams for decay X(4350) — > //^(crossing diagrams should be included). 

Compared to the electromagnetic case, the expression for the matrix element of strong decay is 
more complicated because the constraint from the gauge invariance is released. When X(4350) is 
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regarded as a scalar resonance, the matrix element for the strong decay of X — > V a (qi)Vp{qi) can 
be written as 



G 



9af3 + F : 



Q2 a qip 



£«(<?! )e/?(<?2 



(43) 



x s ^ Vl v 2 so expression 



?1 • <?2 

One can show that, when the gauge invariance is imposed, G Xa _ yV v 
(|43p becomes (|39p . Similarly, without the constraint from the gauge invariance, in the tensor case, 
one can write the matrix element for the strong decay of X^ v — > V a (qi)Vp(q2) as 



n(2) 



(3) 



+F x l l^ ViV2 (g^q u q2 a + gvpq^a) + F^l^ ViV2 q^q u q2 a qi(} e^ v (p)e a (q 1 )e/3(q 2 ) . (44) 



h(5) 



One can prove that when the final vector mesons are both massless particles and the gauge invari- 
ance is imposed the following relations can be reduced 



7 (3) 

X T -*ViV 3 



1 



p(2) 



F 



(5) 



1 ^(2) 



X T -+ViVi 



^X T ->VlV 2 2 (gi • 92)2 X t^ViV 2 ' 



91 • q2 

So expression (|40p for electromagnetic decay matrix element can be yielded. 
With the help of (|42p one can get the analytic forms for the strong decay as 

1 



(45) 



r s (x j/^) 



r T (x -> j/m) 



16^777^ 



2 + 



2G 



1-W 



+ F: 



1 

OJ 



80^777^ 



(46) 



i>j=i 



where u = q\ ■ q2l{jn^m^) = (m x — — m|)/(2m^777^). A is the Kallen function and Cjj are 
functions of the relevant masses of initial and final states which will be given in Appendix. 



IV. NUMERICAL RESULTS AND DISCUSSIONS 

With these discussions, the numerical calculation can be performed via standard loop derivation. 
Since the magnitude of A x is unknown, we vary its magnitude from 0.5 GeV to 4.0 GeV to find 
its physical region where the data can be understood. In our estimate, we use the cental value of 
the total width, i.e., Tx = 13.3 MeV. And, because it is difficult to determine the relative signs 
between and C^ v and other terms, we will consider two cases when we do our numerical 
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TABLE I: Our numerical results in case of positive sign of Eq. (|38"|) . 





A x (GeV) 


Branch product (eV) 


r str (KeV) 


r em (KeV) 


0++ 
2++ 


0.5 - 0.7 
1.1 - 1.8 


2.19- 10.26 
1.24-2.28 


100.9 - 174.5 
285.3 - 973.5 


0.29 - 0.78 
0.03-0.09 


TABLE II: Our numerical results in case 


of negitive sign - 


of Eq. ([23). 


J PC 


A x (GeV) 


Branch product (eV) 


r str (KeV) 


r em (KeV) 


0++ 

2++ 


0.5 - 0.6 
1.0 ~ 1.9 


7.21 - 12.74 
0.66 - 2.42 


373.6 - 391.0 
166.0- 915.1 


0.26 - 0.43 
0.02-0.19 



calculation, i.e., the last two terms of Eq. (j38H give positive and negative contributions to the total 
Lagrangian. Our results are summarized in Tables. [J and [Til 

From the numerical results, one can see that the possibility that A" (4350) is a molecular state 
of D* Q D* can not be ruled out in our model. In the case that AT (4350) has quantum numbers 
J PC = ++ , the physical region of A x is smaller than the tensor resonance case which means the 
size of scalar A (4350) is bigger than the tensor one. 

We would like to point out that, because we used the minimal values of the theoretical calculation 
of coupling constants g D , and g D , our final results about the partial widths can be regarded 
as lower limit. This is an ambiguity of the present calculation. In fact, the best way to determine 
these coupling constants is from data, but because of the precision of the data, we cannot along this 
way. When the magnitudes of coupling constants g D , and g D iDsi are improved, the theoretical 
results of the product of the two-photon decay width and branch fraction to J/ip(f) should be larger 
than the present conclusion. In this case, compared to the tensor AT (4350), the typical region of 
A x for scalar resonance can be reduced to an unphysically small region so one can first rule out 
the possibility of a scalar molecule. 

Another ambiguity in our calculation of the product of the two-photon decay width and branch 
fraction to J/ifj(j> is from the total width of AT (4350). Here we apply the central value, i.e., Tx = 
13.3 MeV. When a larger total width is applied, the physical region of A x can be enlarged. But 
this does not effect the partial widths for strong and electromagnetic decays we predict in the 
corresponding region of A x . 

Finally, we conclude that, with the present data and in the framework our model, A"(4350) can 
be interpreted as D* D* molecule. 
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Appendix A: Explicit forms for the Functions Cy 

In this appendix, I will present the coefficients Cjj in formula 

5A 2 (A + 12m?m|) 

A(^5A(m 2 + m 2 + m|) + 2Am 2 m\m\ 
5A 2 ((m 2 -m 2 ,) 2 
5A 2 ((m 2 -m 2 ) 2 
5A 3 (m 2 — m\ — m|) 
5A 2 + 44m 2 (m 2 + m|)A + 528m 2 m 2 .m 4 
A(5A + 44m 2 m 2 )(m 2 — m 2 + m 2 ) 
A(5A + 44m 2 m|)(m 2 - m| + mf) 
5A 2 (m 4 - (m 2 - m 2 ,) 2 
A 2 (5A + 44m 2 m 2 
5A 2 ^m 4 — (m 2 — m 2 ,) 2 
5A 3 (m 2 - m| + m 2 ) 
A 2 (5A + 44m 2 m 2 , 
5A 3 (m 2 - m? + m|) 
5A 4 







1 p 




24m 4 m 2 m| - 






-1 f 




12m 4 m 2 m 2 . 


W3 




1 [ 




12m 4 m 2 m2 - 


Ui4 




-1 p 




12m 4 m 2 m2 - 






1 |- 




2Am A m\m\ - 


L-22 




1 |- 




2Am A m\m\ - 


<^23 




-1 |- 




12m 4 m 2 m2 - 


<-^24 




1 p 




12m 4 m 2 m 2 - 


L'25 




-1 p 




2Am A m\m\ - 


^33 




1 [ 




2Am A m\m\ - 


C34 




-1 |- 




12m 4 m 2 m2 - 


C35 




1 p 




2Am A m\m\ . 


C44 




1 [■ 




2Am A m\m\ - 


C45 




-1 |- 




24m 4 m 2 m2 - 


C55 




1 [■ 




96m 4 m 2 m2 - 



(Al) 



where A = A(m 2 , m 2 , m|) is the Kallen function and m = m x , mi = m^,, mj 



m w 
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